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Abstract 

An effective spin concept is introduced to examine the mathematical and physical analogy 
between phase coherent charge transport in mesoscopic systems and quantum operations on 
spin based qubits. When coupled with the Bloch sphere concept, this isomorphism allows 
formulation of transport problems in a language more familiar to researchers in the field of 
spintronics and quantum computing. We exemplify the synergy between charge tunneling and 
spin qubit unitary operations by recasting well-known problems of tunneling through a delta 
scatterer, a resonant tunneling structure, a superlattice structure, and arrays of elastic scatterers, 
in terms of specific unitary operations (rotations) of a spinor on the Bloch sphere. 



PACS: 73.23.Ad,03.65.Nk,03.67.Lx,76.60.-k 
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I. Introduction 



Two major areas of research in condensed matter physics are phase coherent charge transport 
in mesoscopic structures HI 121 El and spin based quantum computing 0] El [6l Ul El HI HQl ITTTl . 
These two areas are seemingly disparate and until recently have evolved independently. Some 
efforts to stress analogies between the two fields have appeared recently |[T2l [T3l . In this paper, 
we investigate and develop further the isomorphism between these two areas by introducing an 
effective spin concept to describe phase coherent charge transport through two-dimensional arrays 
of elastic scatterers. 

In tunneling problems, the mesoscopic structure through which an electron tunnels, is char- 
acterized by an arbitrary potential barrier. The transmission and reflection amplitudes are usually 
calculated by the so-called "scattering matrix approach" lfT4l [T5l . The scattering matrix relates 
the incoming (a + ,b~) to outgoing wave amplitudes (b + ,a~) on both sides of a scattering region 
(mesoscopic structure), as shown in Figured! such that 
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where S is the scattering matrix. 

For single-mode transport, assuming an electron incident from the left, 
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whereas, for an electron incident from the right, we have 
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(4) 
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and 



r' 



WiOUT)) = 



(5) 



The tunneling problem is completely characterized by the amplitudes (t, r) or (V, f) depend- 
ing on the direction of incidence of the incoming electron. 

Without any loss of generality, we can always think of the two-component column vector 
\if/{OUT)) as a spinor, since it is normalized in the case of coherent transport. The normalization 
follows from the unitarity of the scattering matrix, i.e., 5 1 5 = /. Furthermore, the spinor \i//(OUT)) 
can be thought as the output of a one-qubit quantum gate whose input is the spinor \if/(IN)) = (1 , 0) t 
or (0, l^where | stands for Hermitian conjugate) depending on the direction of propagation of the 
incident electron. The 2x2 unitary matrix linking the spinors \if/(IN)) and \if/{OUT)) can therefore 
be viewed as the matrix characterizing rotation of a qubit whose initial state was \if/(IN)) and whose 
final state becomes \if/(OUT)). This matrix is also the scattering matrix describing the tunneling 
problem. Herein lies the analogy between quantum logic operation on a spin qubit and coherent 
charge transport in a mesoscopic structure. This paper explores this analogy for single channel 
charge tunneling through a single delta- scatterer, a resonant tunneling structure, a periodic array of 
delta scatterers, and one-dimensional arrays of randomly distributed elastic scatterers. 



Consider the tunneling problem of an electron incident from the left on an arbitrary one- 
dimensional conduction band energy profile E c (x). We refer to the (2x1) column vector \i// l (OUT)) 
in Equation (3) as the effective spin whose components characterize completely the scattering am- 
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plitudes of the tunneling electron. For an arbitrary potential energy profile E c (x), the amplitude 
\ij/ l (OUT)) can be found by successively cascading scattering matrices associated with "subsec- 
tions" within each of which E c (x) is approximated by constant values E cl , E c2 , E c3 ... E cn lfT4l[T5i 
The evolution of the pure state \ip*(OUT)) after crossing a number of subsections can be repre- 
sented using the Bloch sphere concept in which the spinor is parameterized as follows [fT6lfT71 



W(OUT)) = e' y 



cos||0> + sin|e^|l> 



(6) 



where y is an arbitrary phase factor and the angles ((p, 6) are the azimuthal and polar angles, as 
shown in Figure |2] 

In Equation (6), |0) and |1) are the (2x1) column vectors (l,0) t and (0, l) t respectively, 
associated with the north and south poles of the Bloch sphere. They are mutually orthogonal, i.e., 
their inner product < 0|1 > = lfl6l . 

To complete the effective spin picture, we consider the following 2x2 matrix |[T8l 



p = W{OUT)){ifj\OUT)\ = 
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Using this density matrix and the Pauli spin matrices (cr v , cr y , cr z ), the effective "spin com- 
ponents" associated with the spinor \iJ/ l (OUT)) are given by 



(S x ) = ^Tr (pcr x ) = ^ (tr* + rf) = Me (rf ) = Me (ft) , 



(8) 



(S y ) = ^Tr(pcr y ) = ?i(tr* - rf ) = RIm(rf ) = -Elm (ft). 



(9) 
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and 



(S z ) = ^Tr (pcr z ) = \ (|tf - |r| 2 ) = \ (l - 2|r| 2 ) = \ ( 2 |t| 2 - l) . 



(10) 



For an electron incident from the right, |</r r (7A0) = ]!)> an d the density matrix 
p' (=\ifr r (OUT))(ifr r (OUT)\) is such that p' = 1 - p, where p is given by Equation (7) and the com- 
ponents (S x ), (S y ) and (S z ) are just the negative of the values in Equations (8-10). Therefore the 
two spinors corresponding to \if/(OUT)) and \ij/ r (OUT)) are mirror images of each other, corre- 
sponding to a reflection through the origin of the Bloch sphere. This means that \ij/ l {OUT)) and 
\ij/ r (OUT)) are orthogonal, which they must be because the scattering matrix is unitary. 

The unitarity of the scattering matrix also leads to: 



Equation (11) shows that the projection of the spinor in the equatorial plane of the Bloch 
sphere reaches a maximum when \t\ = \r\ = 1/V2. Actually, {S x ) 2 + (S y ) 2 is proportional 
to \t\ 2 (l - |/j 2 ), i.e., the low frequency shot noise power for the tunneling electron |fT9l . Since 
((S x ), (S y ), (S\,)) are proportional to the components of the spinor \if/(OUT)) on the Bloch sphere, 
Equation (12) simply states that the spinor stays on the Bloch sphere during cascading of scattering 
matrices. This is expected for the case of coherent transport. The angles (y, 9, <p) appearing in the 
generic expression of the spinor (or qubit) in Equation (6) can be expressed in terms of the phases 



(S x ) 2 + (s y ) 2 = n 2 \t\ 2 (i-\t\ 2 ), 



(ii) 



and 



(S x ) 2 + (S y ) 2 + (S z ) 2 = h 2 /4. 



(12) 
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and magnitudes of the reflection and transmission coefficients: 



\t\ 

\ r \ e K<t>R-<pT) 

where 4> R and tp T are the phases of the reflection and transmission amplitudes, respectively. 
We get 
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\t\ = cos -, 



WiQUT)) 



t 




\t\e i<t,T 












r 




\ r \ e i<p R 





e 



|r| = sin - = Vl - W, 



and therefore, 
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2 ' \\t\ 
Equations (8-10) are therefore equivalent to 



n 




sin 8 sin ip, 



(20) 



and 




cos 9, 



(21) 



Equations (8) and (9) clearly show that the averages (S x ) and (S y ) contain more information 
than the sample conductance alone. The latter depends only on the magnitude of transmission \t\ or 
reflection \r\ in the Landauer picture [20], whereas (S x ) and (S y ) depend on the phase relationship 
between t and r as well. The phase relationship is a strong function of the energy of the incident 
electron. At non zero temperature, there will be a thermal spread in the energy of the incident 
electron which will lead to a rapid wash out with temperature of the components (S x ) and (S y ), 
i.e., the off-diagonal components of the density matrix p. Note that while (S x ) and (S y ) depend 
on the off-diagonal components of the density matrix and are very energy sensitive, (S z ) depends 
only on the diagonal components of the density matrix and is much less energy sensitive. 

II. 1 Quantum computing gate analog 

The 2x2 unitary matrix or quantum computing gate U QG which relates \ifr(OUT)) and \if/(IN)) 
on the Bloch sphere has the general form [fT6l 



where (a,/3, rj, are real numbers and the R y and R z are the 2x2 matrices associated with rotations 
of the spinor on the Bloch sphere about they and Taxis, respectively. Using the fact that R y (rf) = 



U QG { a ,p,rj,0 = e ia R z (j3)Ry(ri)Rz(a 



(22) 
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e _, s°> and/? z (£) = g-'t^llT6ll. we obtain: 



U QG (a,j3, 



For |^(/A0) = |0>, we have 
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\i//(OUT)) = U QG (a,fi,rj,O\0) 



e <-H)si„l 



(24) 



which is the special case of a spinor on the Bloch sphere in Equation (6), corresponding to 



a = y = (f) T , 
i] = 6 = 2tan~ l 
fi = -C = <p = (p R - <p T - 



(25) 



Hence, from a quantum computing perspective, the analytical expression for Uq G is identical 
to the scattering matrix used to described the tunneling problem and is given explicitly by 



UqgWt, 6, \t\) = e'^ T R z ((f>R - (p T )R y \2tan- 1 



R z ((p T - cp R ), 



(26) 



This last equation helps visualizing coherent charge transport (or tunneling) through specific meso- 
scopic devices as a successive set of rotations of the effective spin on the Bloch sphere, as will be 
illustrated in the numerical examples in section III. 

In the next section, we provide several examples to illustrate the effective spin concept. 



III. Examples 



IIL1 Scattering across a single delta scatterer 
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We first determine the quantum computing gate analog of a simple delta scatterer of strength 
Vi6(x) for which the reflection and transmission amplitudes are easily shown to be 



t' = t= —— = , (27) 

ik - ko ik - 1 



and 



r = r = = — , (28) 

ik - k {) ik - l 



with k = k/k , ko = m*Vi/ti 2 and k = V2 ^* g , where E is the kinetic energy of the electron and 
is its effective mass. 

The magnitude and phase of t and r are therefore 



m* 



\t\ = -^=,<Pr = -tan-Mi], (29) 
ik 2 + 1 w 



and 



\r\ = 1 , <p R = tan" 1 (^) - n. (30) 
W + 1 

The spinor \i// l (OUT)) for this simple problem is given by Equation (6), where 



n 

(p = (p R - (f> T = --, (31) 



and 



6» = 2tan" 1 (l/I). (32) 
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The equivalent quantum computing gate is characterized by unitary matrix U QG given by 



U QG = e ^R z (-?-)R y (6)R z (^) = e^R x (-9), 



(33) 



where R x is the matrix for spinor rotation around the x-axis lfT6ll . For low incident energy, 8 = n 
and it monotonically goes to as the energy of the incident electron increases. According to 
Eqns. (19-21), the spinor \if/\OUT)) sweeps only a very limited portion of the Bloch sphere, i.e., 
the semi-circle in the y-z plane, going from the south to north poles clockwise as the energy of the 
incident electron increases. The spin components of \i]/ l (OUT)) along the x, y, and z axes are given 

by 



(S x ) = 0, 



(34) 



(s y ) = 4 



( 2k 
k 2 + l 



(35) 



and 



2U 2 + 1 



(36) 



For instance, when k = 1, \ifr'(OUT)) is in the equatorial plane of the Bloch sphere, along the 
y-axis. In this case, 9 = |, and the matrix Uq G is given by 



U QG = e-*R z (--)R y (-)R z (-) = e~S {-^)cr/HS (-), 



(37) 



where 



S(S) 



1 

e iS 



(38) 
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and 



V2 



1 1 
1 -1 



(39) 



are the general phase shift and the Hadamard matrix, respectively, extensively used in the theory 
of quantum computing [T6l . 

IIL2 Scattering through a delta-scatterer in a region of length a 

Next, we consider the scattering problem across a region of length a containing a delta scat- 
terer at location xq. The corresponding scattering matrix can be easily derived. The location of the 
spinor \ifj l (OUT)) on the Bloch sphere is described by azimuthal angle 6 given in Equation (21) 
and polar angle 



-7T 



<p = — -k(a-2x Q ). 



(40) 



The average values of the effective spin components are given by 



(S x ) = 

(Sy) = 

(S z ) = 




sinfc(2xo - a), 
cos k (2xq - a), 



\k 2 + l) 



(41) 
(42) 
(43) 



In this case, (S x ) is non-zero unless xq = |, i.e., unless the potential energy profile in the 
device is spatially symmetric. For a fixed value of the incident wavevector, the spinor \if/(OUT)) 
moves on a circle parallel to the (x, y) plane. If a is selected such that ka = n, (p increases linearly 
from -y to | as xq varies from to a, i.e., the Bloch vector associated with the spinor sweeps 
the entire plane defined by the component (S z ). According to Eqns.(41) and (42), if ka = n, the 
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average value of (S x ) and (S y ) are equal to zero when we average over the impurity location x . 
This is an important ingredient in the theory of localization in ID arrays of scatterers, as will be 
discussed later. The quantum computing gate U QG analog of this tunneling problem is given by 

U QG = e^R z (-| + k(2x Q - a))R y (6)R z (k(a - 2x Q ) + ~) . (44) 

Since 6 is still given by Equation (21), a spin flip from the south to north pole is only possible 
if we increase the energy of the incident electron to infinity. The energy cost for the spin flip is 
drastically reduced if we have two or more delta scatterers, as discussed next. 

IIL3 Scattering across a resonant tunneling structure 

We consider the scattering problem across a resonant tunneling structure consisting of two 
delta scatterers of equal strength Vj separated by a distance a. In our numerical simulations, we 
use V t = 0.3eVA and a = 50 A. Figure [3] is a plot of the transmission coefficient T as a function 
of the reduced wavevector k. The first two resonances (at which T = 1) occur at k « 12.5 and 36. 
The corresponding variation of the phase angles (<p, 6) for the spinor \ifz\OUT)) are displayed in 
Figure SI The angle 6 reaches its minimum value of zero at the resonances when there is a sudden 
jump in <p. When viewed as a quantum computing gate, an RTD is more efficient when operated 
over the range Ak indicated in Figure H] since it allows a full swing in 6 from to n, whereas the 
swing in 9 is much smaller between the first two and higher resonances. The quantity Ak is much 
smaller than the infinite change in k needed for a single delta-scatterer to realize an inverter, as 
discussed in the previous section. Since T = R for k = k* , 8 = | which is enough to implement the 
Hadamard gate using an RTD. 

The results above can be extended to the case of a superlattice, modeled as a sequence of 
evenly spaced identical delta scatterers. In that case, each resonant state present in the smaller unit 
with two scatterers (RTD) leads to a passband for the infinitely periodic structure. In Figure [51 
we plot the transmission coefficient for a structure consisting of 5 delta scatterers with the same 
parameters as for the RTD described above and with the same spacing of 50 A between each 
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scatterer. The transmission coefficient reaches unity at four vaiues of k in the intervai [5 - 25], 
which is a weif known resuft for finite repeated structures ll23l |24| . Furthermore, the range Ak 
needed to reach the condition T = R is reduced compared to the case of a RTD. As the number 
of periods in the superlattice increase, Ak actually converges to a limit corresponding to the lower 
edge of the pass band of the infinite superlattice. As shown in Figure [6l the angle allows a full 
swing from north to south poles on the Bloch sphere over a range Ak smaller than what is necessary 
for the case of the RTD, and the phase angle <p toggles back and forth between -| and | each time 
a resonance is crossed. 

Figure[5]also shows a plot of the transmission coefficient (curves labeled 1 and 2) versus k for 
two imperfect structures, in which the locations of the five delta scatterers are selected randomly 
and uniformly over each interval of length a. The transmission coefficient is fairly sensitive to k in 
the range of k where the lower pass band will develop for the infinite superlattice. However, the 
transmission curve is fairly insensitive to the imperfections in the superlattice in the same range of 
k. As shown in Figure |6j the angle 6 is also fairly insensitive to imperfections in the superlattice 
but the phase (f is not. The latter result is a compounded effect of multiple reflections between 
impurities and the sensitivity of (f to the exact impurity location in each section of length a, as 
discussed earlier. 

III. 4 Scattering through a periodic array of delta scatterers 

The scattering matrix elements for ID periodic system (or superlattice) can be calculated 
exactly B25l|26l|271. The transmission amplitude is found to be 



tN — 



D N 



(45) 



and the reflection amplitude is given by 



r N = -i 



k e i(N - l)ka sin(A^a) 

i , 

k D N sm(fia) 



(46) 
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where 



D N = e 



iNka 



cos(N]3a) + Urn 



sin(Nf3a) 



(47) 



k J sin(/fa) 

and k = \ y/2m*E, k = m*Vi/fi 2 , a is the distance between adjacent scatterers, and fi is the quasi 
momentum. It is the solution of the transcendental equation: 



ko 

cos(/3a) = cos(fca) + —sin(ka). 

K 



(48) 



Using Eqs.(8-10), it can be shown that < S x >= and 



<S y > = -hRe(rt) = - 



sin(Nf3a) k 



i , / ^0^2 sin (Nfia) sin(fla) k 

y k ' sin 2 (J3a) 



and 



<Se> = ?(2| t P-l) = 5 



L ^^k> S in 2 03a) 



- 1 



(49) 



(50) 



In the case ofN= 1, we get back Eqs. (35) and (36) of section DLL In the case of delta 
scatterers, incident energies for which 



sm(N/3a) 
sin(J3a) 



= 0, 



(51) 



correspond to points of unity transmission which occur at values of the quasi-momemtum in the 
first Brillouin zone 

7777 

(52) 



nn 

B n a = — 



with (/i = \,...N- 1). 

At these values, < S x >=< S y >= and < S z >= \. 

III.5 Transport through random arrays of delta scatterers 

The analysis of the previous section was extended to a large number of delta scatterers of strength 
Vi6 - (x l + (i - l)a)), where Vj is selected to be 0.3 eVA and x l Q is the location of the i th impurity 
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located in the interval [(?' - \)a, ia] Each impurity location is generated using a uniform random 
number in each interval. The length of each subsection is set equal to 237 A and the wavevector of 
the incident electron k = ^ 2l f E , is selected such that ka = n, for an incident energy £ of 10 meV 
and m* = 0.067m , the electron effective mass in GaAs. 

Figure |7j is a plot of the phase angle of the spinor \ij/ ! {OUT)) versus the number of sub- 
sections (TV) crossed. The two top curves are 6 versus /V for two specific impurity configurations. 
The curves show regions where 6 decreases as /V increases which corresponds to an increase in 
the conductance of the array. This decrease in 6 as /V increases is quite pronounced for one of the 
two impurity configurations, for /V < 20. A plot of the average value of 6 over an ensemble of 
10 5 samples is shown as the curve labeled 8 in Figure [7] The quantity = n/2 for /V « 23. This 
corresponds to a conductance of e 2 /h, as shown in Figure[8l and to an elastic mean free path equal 
to 23 x 237A ~ 0.55/im. 
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IV. Conclusions 



The effective spin concept examined in this paper offers an alternative description of phase 
coherent charge transport through mesoscopic systems in terms familiar to researchers in the field 
of spintronics and quantum computing. As illustrated in this paper, the effective spin formalism 
provides a pedagogical approach to simple scattering problems and also to the phenomenon of 
localization in random arrays of elastic scatterers. 

In the past, the effective spin concept has been used to describe the spatial correlations be- 
tween reflection and transmission amplitudes of polarized photon beams from a combination of 
beam splitters, mirrors, and interferometers [T3T1I321I331I . More recently, the effective spin concept 
has been used to examine the critical problem of entanglement between channels associated with 
propagating modes in mesoscopic systems, as reported in recent experiments by Neder et al. lfT2ll 
and their theoretical interpretation by Samuelson et al. [|34l . 
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Figure Captions 



Figure 1: The tunneling problem and its quantum computing gate equivalent. The scattering ma- 
trix associated to a device relates the incoming (a + , b~) to the outgoing (a~, b + ) wave amplitudes. 
It can be viewed as the matrix representing the rotation of a qubit from the initial state \iJ/(IN)) to 
the final state \ifr(OUT)). 

Figure 2: Bloch sphere representation of the effective spin (qubit) \if/{OUT)). The radius of the 
sphere is equal to 1 . 

Figure 3: Transmission (T = \t\ 2 ) and reflection (R = |?| 2 ) coefficients versus reduced wavevector k 
of electron incident on a single delta scatterer. The expressions for \t\ and \r\ are given by Equation 
(51) and (52), respectively. 

Figure 4: Plot of the phase angles {tp, 0) associated to the spinor \ifj(OUT)) describing tunneling 
through a resonant tunneling structure as a function of the reduced wavevector k = k/ko, where 
k = \ y/2m*E, E is the kinetic energy of the incident electron in the contact; k = m*Vil% 2 , and 
Vj is the strength of the delta scatterer. The two delta scatterers are separated by 50A and have a 
strength Vj = 0.3eVA. Ak is the minimum wavevector (in reduced units) needed to realize a spin 
flip from the south to north poles on the Bloch sphere. The zeroes in 6 are the locations of the 
energy resonances. 

Figure 5: Transmission (T) and reflection (R) coefficients versus reduced wavevector k of the 
incident electron for a superlattice modeled as five delta scatterers of Vj separated by a distance 
a (Vi = 0.3eVA and a = 50 A). The curves labeled "1" and "2" are T versus k for two imperfect 
superlattices, i.e., for two arrays of 5 delta scatterers whose positions are selected uniformly over 
each interval of length 50A. 

Figure 6: Reduced wavevector dependence of the phase angle (</?, 9) associated to the spinor 
\i//(OUT)) describing tunneling across an array of five delta scatterers separated by 50A and with 
a scattering strength Vj = 0.3eVA. The zeroes in 6 are where the transmission through the super- 
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lattice reaches unity. Also shown as dashed lines are the angles (<p, 6) through two random arrays 
of elastic scatterers. 

Figure 7: Evolution of the angle 6 for the spinor \i]/{OUT)) on the Bloch sphere as a function of 
sample length for two different arrays of elastic scatterers (two top curves). The smoother curve 
represents the average of 6 calculated over an average of 10 5 arrays with the locations of each 
individual scatterer varied uniformly across each subsection of the array. The elastic mean free 
path A e i (in units of subsections crossed) is where 6 = | . 

Figure 8: Plot of the average over an ensemble of 10 5 impurity configurations of the conductance 
as a function of the number of impurities crossed in the sample. Also, shown is the value of the 
classical conductance calculated neglecting the effects of multiple reflections between scatterers. 
The elastic mean free path A d (in units of subsections crossed) is where the Landauer conductance 
reaches a value of e 2 /h. 
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Figure 6: J. Wan et al 
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